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Abstract. High quality observations of the low-degree acoustic modes (p-modes) exist for almost two complete 
solar cycles using the solar spectrophotometer Mark-I, located at the Observatorio del Teide (Tenerife, Spain) 
and operating now as part of the Birmingham Solar Oscillations Network (BiSON). We have performed a Fourier 
analysis of 30 calibrated time-series of one year duration covering a total period of 15 years between 1984 and 
1999. Applying different techniques to the resulting power spectra, we study the signature of the solar activity 
changes on the low-degree p-modes. We show that the variation of the central frequencies and the total velocity 
power (TVP) changes. A new method of simultaneous fit is developed and a special effort has been made to 
study the frequency-dependence of the frequency shift. The results confirm a variation of the central frequencies 
of acoustic modes of about 0.45 fiRz, peak-to-peak, on average for low degree modes between 2.5 and 3.7 mHz. 
The TVP is anti-correlated with the common activity indices with a decrease of about 20% between the minimum 
and the maximum of solar cycle 22. The results are compared with those obtained for intermediate degrees, using 
the LOWL data. The frequency shift is found to increase with the degree with a weak ^-dependence similar to 
that of the inverse mode mass. This verifies earlier suggestions that near surface effects are predominant. 
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1. Introduction 

Understanding the observed solar variability is one of the 
major goals of solar physics. Because the frequency shifts 
of solar p-modes are known to be very sensitive to the solar 
activity cycle, the analysis of helioseismic data has been 
used to track those physical processes which underly the 
origin of the cyclic changes observed at the solar surface. 
Helioseismology based on low-degree p-modes is necessary 
to look for potential structure or dynamic changes in the 
deep interior. 

The first report of frequency shifts of the low -degree 
p-modes was given by Woodard fc Noyes (1985 ). Using 
ACRIM data, they found that the few observed £=0 and 
1 modes presented a change in the central frequency of 
0.42±0.14 /xHz in average during the declining phase of cy- 



the south pole, and by Palle et al. (1989) using a long set 
of data from the Mark-I instrument at Observatorio del 
Teide covering the full cycle 21 (1977-1988). 



Subsequently, [Rcgulo ct al. (1994[ ) used Doppler ob- 
servations collected from the maximum of cycle 21 to the 
falling phase of cycle 22 (1980-1993), obtained with Mark- 
I instrument, to calculate monthly frequency shifts. They 
showed that, for all low degree acoustic modes, there is 
an important frequency shift of 0.52±0.02 /iHz correlated 
with solar activity. In addition, the amplitude of these 
variations is different when £—1,3 and £=0,2 are consid- 
ered separately. The odd modes present, on average, a 
change of 0.58±0.06 /LtHz; the even ones show a full shift 
of only 0.33±0.06 /^Hz. 

Two other important properties of the low-degree 



cle 21 (1980-1984) . These results were confirmed by |Fossat| changes have also been pointed out recently [Anguera 



et al. [i^»Y|j by comparison witn observations made at q^^^^ ^1. (1992]) observed the frequency dependence 
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of the freque ncy shifts for the low -degree p-modes, (later 
confirmed by Chaplin ct al. 1998) in agr eement with the 
earlier resul ts of Woodard et al. (1991 ) for intermedi- 
ate degrees. Jimenez-Reyes et al. (1998 ) found that the 
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frequency shifts, when plotted against an activity index, 
show a hysteresis behavior rather than a simple linear cor- 
relation. This result was interpreted as part of structural 
changes associated with the solar activity which are tak- 
ing place in the Sun. This must be confirmed by more 
observations but the interpretation of these later results 
as being partly due to structural changes in the interior 
associated with the solar activity has been found to be 



a complex problem. Recently, Moreno-Insertis & Solanki 
(2000|)]have studied in detail the signature left on the low 
degree p-mode frequencies by the surface solar magnetic 
activity. Whether these changes are taking place only close 
to the surface or not is not completely clear and one of the 
requirements to address this question is to get precise and 
reliable measures of the low degree mode parameters for 
a long period of time. 

The spectrophotomer Mark-I, has been collecting so- 
lar observations for almost two complete solar cycles. The 
available database for low-degree p-modes, probably the 
longest in duration and the most stable, is used in the 
present work to analyze the signature of the solar cycle 
in the mode parameters and to parameterize the observed 
frequency shifts as a function of various classical solar in- 
dices. 

In the following section, the essential steps of the data 
reduction leading to the yearly spectra are presented. In 
general, the frequency shift between a time U and a time 
to taken as reference, can be written as a function of the 
frequency and the degree, i.e. Siy{t — to,v,i). In the fol- 
lowing the reference time is 1986 which corresponds to a 
minimum of solar activity. In Sect. ^, the integrated fre- 
quency shift Ai^' =< 5v{ti — to, ly^i) > is analyzed where 
the brackett indicates an average for all observed low de- 
grees and frequencies between 2.5 and 3.7 mHz. Two tech- 
niques are proposed to measure the frequency shift from 
the cross-correlation function between power spectra of 
time-series created at different solar activity level. In ad- 
dition to the frequency shift, the second technique allows 
us to study the time variation of the total velocity power 
(TVP) which are presented in Sect. ^. Then, in Sect. ||, 
we focus on the study of the frequency dependence of the 
frequency shifts. Again, two different techniques are used. 
The first one consist in simply cutting the spectra in band 
of 135 /iHz before computing the cross-correlation func- 
tions. The second one is a new procedure developed here 
and called simultaneous fitting: all the yearly spectrum 
are fitted at the same time assuming that the time depen- 
dence of the mode frequencies can be described as a linear 
function of the radio flux at 10.7 cm FIq, taken as solar 
activity index i.e.: 

Si^it, - to, y, €) = ^v{v, £){Flo ~ F°o), 



(1) 



where F°q represents the radio flux at the 1986 solar min- 
imum. In addition, we define 5i'{v) as the frequency shift 
per radio flux unit averaged over L We note that if we 
had used shorter time-series and a magnetic index in- 
stead of the radio fiux to parameterize the time depen- 
dence of the frequency shift, a more complicated formu- 



lation would probably have been needed in order to take 
into account the hysteresis behaviour found when mag- 
netic indices are plotted versus the frequency shift during 
the cycle (Jimenez-Reyes et al. 1998). 

In order to check the ^-dependence of the frequency 
shift, we use the fact that pairs of low-degree p-modes 
with the same parity have "almost" the same frequency 
and, that they are equally spaced in frequency. This allows 
us to provide in Sects. ||, ^, ^ not only the ^-averaged 
quantities defined above (Ai/, TVP, 5v{v)) but also the 
quantities related to even and odd modes separately i.e. 
respectively: Aiyo,2, TVPo,2, 5v{v)q^2 and Ai^i,3, TVPi,3, 
8ii{i>)i^^. Finally, in Sect. ^ the ^-dependence analysis of 
the frequency shift is completed by comparing the results 
with those obtained at higher degrees [1=1, 99) using the 
LOWL database. 



2. Observations and data analysis 

The data used in this work come from the observations 
carried out at the Observatorio del Teide between 1984 
and 1999. The observations consist of daily measurements 
of the solar radial velocity obtained with the Mark-I res- 
onant scattering spectrophotometer. This instrument has 
been sited at the Observatorio del Teide since 1975. After 
some hardware updates in 1984, the experiment has been 
running without interruption other than bad weather and 
instrumental failures. The data reduction process is ex- 



plained in more detail elsewhere (van der Raay et al. 1985 



Palle et al. 1986, 1993). Briefly, the data are corrected from 



the annual scan (Earth's orbit around the Sun) of the non- 
linear solar line shape, and calibrated by fltting the known 
daily velocity of the observatory; only two parameters are 
fitted, taken to be the same over the all 15 years analyzed 
here. Then, the daily residuals are joined in consecutive 
360 days leading to a total of 30 time-series with 6 months 
in common between consecutive series. Although all series 
are not independent, they show very similar duty cycles 
(around 25%). Finally, the corresponding power spectra 
were calculated for every time-series using a traditional 
Fourier analysis. All the spectra show the peaks of the low- 
degree p-modes with f < 3. The typical sideband structure 
appears in the spectra at kjj— 11.57 fiHz as a direct con- 
sequence of the observing window achieved from just one 
station. 



3. Frequency integrated frequency shifts 

Once the power spectrum of each time-series has been 
calculated, the integrated frequency shifts A;y* are deter- 
mined. This is done by computing the cross-correlation 
p^{vj) of each power spectrum i with the power spectrum 
of the time-series covering 1986 taken as reference. As we 
show in Sect. 5, the dependence of the frequency shift 
for low and intermediate-degree p-modes appears to be 
essentially null below 2 mHz, whereas at high frequency 
(above 3.7 mHz) the frequency shift is expected to drop 
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quickly. The analysis of the integrated changes is lim- 
ited to the modes between 2.5 and 3.7 mHz. considering 
a larger interval could bias the results by including the 
strong high frequency variations that are less accurately 
determined and have possibly different physical explana- 
tion dGoldreich et al. 1991j ). 

In order to calculate the position of the cross- 
correlation main peak, two different methods have been 



used. The first was used in previous analyses (Regulo et al 
1994; Jimenez- Reyes et al. 1998 ) and takes the maximum 
of a second-order polynomial fitted to the logarithm of the 
cross-correlation function in an interval iter around the 
main peak (where a is the second-order moment). This 
cross-correlation is calculated starting at an appropriate 
lag for which the function is symmetric, and which is ob- 
tained by calculating the third-order moment. This proce- 
dure directly provides a value of the mean frequency shift 
between the £ < 3 p-modes and the corresponding values 
at solar activity minimum, the chosen reference. 

The second method, introduced here, is based on the 
shape of the cross-correlation function. Assuming that 
each oscillation mode can be modeled by a damped har- 
monic oscillator, each peak in the spectra has a Lorentzian 
profile and the correlation function has also a Lorentzian 
profile. In order to improve the determination of the fre- 
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Fig. 1. Differences between the integrated frequency shifts 
as inferred from the cross-correlation functions by the two 
methods explained in the text i.e. (1) fit by a polynomial 
close to the maximum and (2) fit by a Lorentzian profile 
including side bands. The differences remain small and not 
significant. 

modeling of the p-mode excitation and damping mecha- 
nisms. It uses only information contained in the spectrum 
while the second method assumes symmetric Lorentzian 
profiles. There is some evidence that the shape of the 
peaks in the power spectra are slightly asymmetric (Thiery 



quency shift, the model also take into account the presence et al. 2000). In order to check the efficiency of both meth 



of sidebands with amplitude located at ±11.57 /iHz. 
Thus, our model of the cross-correlation function between 
the spectrum i and the reference spectrum can be written 
as: 



k=-\ 



A*(rV2)2 



B\ (2) 



where the parameters to be fitted are: 



— A'', the amplitude of the central peak; 

— F', the linewidth of the Lorentzian profile; 

— A^*, the average frequency shift in the chosen fre- 
quency interval; 

— 5% the constant background level; 

— , the ratio of the sidebands to the central peak {(3—1 
for fc=0 ); 

— ko = 11.57/xIIz is the constant separation of the side- 
bands and it is the only fixed parameter in the fit. 

The best estimation, in the least square sense, of the pa- 
rameters related to each time-series i, is obtained by min- 
imizing the following quantities: 



N 



(3) 



where N is the total number of frequency bins in an in- 
terval of ±20 /iHz around the main peak of the cross- 
correlation function. We applied a Levenberg-Marquard 
method (Press et al. 1992) but any other minimization 
routine may be used. 

We note that the first technique is more objective in 
the sense that it does not require a physical hypothesis or 



ods and the influence of the duty cycle in the final result, 
we analyzed the frequency shifts in periods of 36 days for 
which the duty cycle varies. The differences between the 
two methods remain in general within the error bars. The 
only significant differences are found for time-series with 
very small duty cycle (around 10%) but, for the yearly 
time-series analysed here, the duty cycle is moderately 
high and quite stable (around 25%) from year to year. 
Regarding to the observed asymmetry of the p-modes, 
they are not thought to be important for this analysis. 
The second method takes into account the known distance 
from the sidebands to the main peak (±11.57 /iHz) which 
improve the results. Moreover this method allows to study 
not only the frequency shift but also the TVP by providing 
amplitudes and widths of the cross-correlation functions 
(see Sect. Therefore hereafter only the second method 
will be considered but, for completeness. Fig. |l| shows the 
marginal differences between the frequency shifts obtained 
by applying both methods to the yearly time-series. 

As mentioned in the introduction, the asymptotic the- 
ory predicts that, for low-degree p-modes, pairs of modes 
with alternately odd and even degrees are equally spaced 
in frequency with a separation of about 67 /iHz i.e. half of 
the so-called big separation (e.g. Deubner & Gough 1984). 
The contributions of even (£=0,2) and odd degrees {£=1,3) 
to the integrated frequency shift can therefore be sepa- 
rated simply by applying a mask to the spectra before 
using the procedure explained above. It should also be 
noted that, for full disk observations, the contribution of 
the two modes of a pair is not the same, due to the geom- 
etry of the modes at the surface. In the case of the odd 
modes, the frequency shifts come essentially from £=1 due 



4 



Jimenez- Reyes et al.: Analysis of the solar cycle and core rotation. I 



Table 1. Correlation coefficients between different solar indices and the yearly frequency shift, rp is the Pearson linear 
correlation coefficient, rs the Spearman rank correlation coefficient and Ps is the probability of having no correlation. 

Index Au Avo 2 Aui 3 



rp rs n rp rs rp rs 

"R^I 094 088 9 10""^ 076 0?70 TW^ O90 085 TW^ 

Fio 0.94 0.87 4 10"^° 0.77 0.72 7 10"*^ 0.90 0.83 1 10"* 

KPMI 0.90 0.86 8 10-^° 0.79 0.71 1 10"^ 0.85 0.84 8 lO"'' 

MPSI 0.94 0.88 5 10"^° 0.83 0.76 2 lO"'* 0.88 0.85 1 10"* 

TSI 0.89 0.89 3 10"** 0.78 0.78 1 10"^ 0.83 0.83 2 10"^ 

He 0.94 0.88 9 10"" 0.79 0.74 2 10"*^ 0.89 0.83 2 10"* 
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^ 0.4 
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Fig. 2. Time variation of the integrated frequency shift for 
low-degree p-modes is plotted in the upper panel (black 
dots), where the radio flux at 10.7 cm is also shown (full 
line). In the sub plot, the results for the even and odd 
degrees separately are presented. The averages of these, 
are plotted in the main panel (squares). The lower figure 
shows the good linear correlation between the frequency 
shift and the radio flux. 



to the high ratio in sensitivity between £=3 and 1 (~ 0.1) 
whereas, in the case of the even modes the sensitivity ratio 
is close to one and therefore the contributions of £=0 and 
2 are nearly the same. 

Figure § illustrates the frequency averaged frequency 
shift between 2.5 and 3.7 mHz. The results corresponding 
to ^=0,2 and £=1,3 have been plotted in the sub plot at 



Table 2. Intercept and slope of the frequency shift ex- 
pressed as a linear function of different solar indices. 



Solar 




Intercept a 


Slope b 


Index 




(nHz) 


(nHz per activity unit*) 


Ri 


Af 


18.25±14.77 


2.56±018 




Auo,2 


-48.08±25.94 


1.98±032 




Ai/1,3 


60.62±21.51 


2.91±026 


Fio 


Au 


-139.07±24.17 


2.69±0.19 




Auo,2 


-173.73±41.78 


2.11±032 




Ai^i^s 


-115.88±36.40 


3.03±O28 


KPMI 


Av 


-98.43±27.42 


22.92±2.04 




Auo,2 


-154.75±37.20 


19.08±2.78 




Aui,3 


-63.15±40.42 


25.34±3.02 


MPSI 


Av 


33.79±14.04 


149.06±10.46 




Auo,2 


-39.06±21.10 


124.22±15.59 




Aui^3 


79.81±22.19 


164.63±16.56 


TSI 


Av 


-528.06±51.63 


386.71±37.80 




Al/0,2 


-437.13±65.43 


32O07±47.90 




Ai/1,3 


-585.31±71.58 


428.67±52.40 


He 


Av 


-438.86±42.45 


10.92±0.73 




Afo,2 


-417.59±73.56 


8.87±1.26 




A!^l,3 


-447.34±67.87 


12.21±1.17 



'Units: nHz; nHz/(10 
nHz nHz mA 



J/s/mVHz), 
respectively. 



nHz G"\ nHz 



the top-right corner and their average is shown in the main 
figure. The solid line represents the radio flux averaged 
over the same periods as the time-series used for this work. 

Although there are few departures from the general 
trend which do not agree with the smooth behavior of the 
solar index, the integrated signal concurs very well with 
the radio flux, which represents here the behavior of the 
solar cycle. The average of both contributions (£=0,2 and 
£=l,3) follows well the integrated signal as expected. The 
amplitude, measured as the straight difference from peak- 
to-peak, for all observed modes is 0.45±0.05 /iHz, while in 
the case of the even and odd degree modes are 0.48±0.05 
and 0.55 ± 0.07 /xHz respectively. Although the difference 
in amplitude between even and odd degree modes seems 
to remain, the frequency shift for the even modes is larger 
than that obtained by Regulo et al. (1994). 

Since the time-series created are one year long, the 
time variation of the integrated frequency shift analyzed 
here informs only on long-term changes. A recent analysis 
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covering different time scales from one to seven months 



using 9 years of BiSON data can be found in Chaplin 
et al. (2001). Average values of the following solar activ 



Lean 



ity indices have been computed over the same one year 
periods than the frequency shifts in order to obtain the 
corresponding correlation coefficients: 

— the International Sunspot Number, R/; 

— the integrated radio flux at 10.7 cm, Fio; 

(both obtained from the Solar Geophysical Data)|^; 

— the Kitt Peak magnetic index (KPMI) extracted from 
the Kitt Peak full disk magnetograms ( [Harvey 198^ ); 

— t he Mount Wilson Magnetic Plage Strength Index 
(lUlrich 199l| ), MPSI; 

— the Total Solar Irradiance, TSI ( Frohhch fc 
199|); and 

— the equivalent width of Hel 10830A averaged over the 
whole solar disk using data from Kitt Peak observa- 
tory. 

The Pearson correlation coefficient rp, which is a mea- 
sure of the strength of the linear relationship between two 
indices, and the Spearman rank correlation coefficient, rs 
which provides a measure of the correlation between the 
ranks of two indices during the chosen period, are shown in 
Table |^ for the frequency shift. The correlation between 
the different solar activity indices themselves have been 



investigated in details by Bachmann & White (1994) 



In addition, the probability Ps of having null correla- 
tion between the ranks of any of the solar indices and the 
frequency shift is indicated. We note that the database 
of MPSI and TSI indices, covers only the first 28 time- 
series. The correlation analysis for these two indices was 
therefore made with two fewer points than for the others. 

The integrated signal shows very high correlation with 
the various solar indices, whereas the frequency shift cor- 
responding to £=1,3 and, more significantly, £=0,2 are 
slightly lower. The general trend of the frequency shift 
corresponding to the even and odd degree modes sepa- 
rately is quite clear but, in addition to the difference in 
amplitude already mentioned, the frequency shifts mea- 
sured for the even modes seem to be sensitive to the solar 
cycle later than the odd ones and the resulting difference 
in phase between the two curves is probably at the origin 
of the lower correlation coefficients found for £=0,2. As 
demonstrated by Moreno- Inscrtis fc Solanki (2000| ) differ- 



ent modes would respond differently at different phases of 
the cycle depending on the positions of the activity (i.e. 
sunspots) on the disk. Aside of the long term differences, 
there are also fluctuations at shorter time scales which are 
different for the two data sets. 

Because of the excellent linear correlation coefficients 
found, the frequency shifts were fitted as a linear function 
of the different solar indices / by: 



Aiy' = a + b- r 



(4) 



In Table | we report the intercepts and the slopes obtained 
for all the solar indices considered, and an example is given 
for the radio flux in the lower part of Figure |[ The slopes 
can be compared with the results shown in Table 1 of 



Regulo et al. (1994) also obtained for low-degree p-modes. 



The three solar indices in common to both works present 
similar slopes; the differences are less than 3 times our 
error bars. Notice that, more or less, all activity indices 
used here produce similar values of the correlation coeffi- 
cient with the frequency shifts leading to the conclusion 
than none is much better, for this purpose, than others. 

The numbers als o agree with those obtained recently 
by Jain et al. (2000 ) for intermediate-degrees (i?=20-100) 
using GONG data. However, intermediate degree modes 
are confined closer to the surface and one may expect them 
to be more sensitive to the activity changes and therefore 
the slope to be larger for those modes than for low-degree 
modes. This is indeed what we found analyzing LOWL 
data ( Jimenez- Reyes et al. 2001 ) but, as we shall see, this 
depends also strongly on the range of frequencies consid- 
ered to calculate the average values. 

The interpretation of the different behavior found for 
£ =1,3 and £=0,2 is not straightforward. In order to under- 
stand better the underlying physics, one may instead look 
at the velocity power variations and at the frequency and 
^-dependences. This is considered in the following section. 

4. Total velocity power variations 

The first observations of the variation of the TVP for all 
measured p-modes was reported by Pallc et al. (1990a| Jb|) 



who found an increase of 30 to 40%, anti-correlated with 



the so lar activity cycle. Afterwards, Anguera Gubau et alj 



1992 ) found similar changes using different analysis tech- 



niques and more data of the same type; they interpreted 
these results as a decrease in the efficiency of the excita- 
tion of such modes at solar activity maximum, since ab- 
sorption of mode power by local magnetic structures (see 
e.g. Bogdan et al. 1993) is a small influence and cannot 



^ Available at |http://web. ngdc.noaa.gov/stp/stp. htm] 



explain such a high ratio. 

We have also calculated the TVP in the spectrum 
which is proportional to the area under the main peak 
of the cross-correlation fimction. Once the main peak is 
fitted to a Lorcntzian profile, the TVP is calculated as 
the amplitude times the width (i.e TVP' = A"- ■ F'). We 
show in Fig. ^, the percentile changes of the TVP com- 
pared to its minimum value. Again, the radio flux is shown 
here as an index of the solar activity cycle. From the fig- 
ure, we are able to see that the variation of the TVP be- 
tween minimum and maximum of the solar cycle is around 
20%. Moreover, the lower figure shows clearly the anti- 
correlation with the solar cycle, the TVP decreases when 
activity increases. The separated contributions of the odd 
and even pairs of modes were also studied (see the inset 
box in Fig. |^) as we did in the analysis of the frequency 
shifts. The results are shown in the sub panel on Fig. |[ 
Notice that the variation of the TVP between extreme 
phases of the cycle for even and odd degrees is approxi- 
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Fig. 3. Time variation of the TVP changes in per- 
cent compare to its minimum value (i.e. 100 • (TVP^ — 
TV Pmin) /TV Pmin) for low-dcgree p-modcs integrated be- 
tween 2.5 and 3.7 mHz. In the sub-panel the TVP corre- 
sponding to £=1,3 and £=0,2 are also shown. The radio 
flux at 10.7 cm calculated for the same periods is plotted 
as a full line in the upper plot and its linear correlation 
with the frequency shift is shown in the lower plot. 




5 10 15 

TVP changes (%) 



25 



Fig. 4. TVP increase in percent against the frequency 
shift. The solid line goes from 1984 to 1989 and the dashed 
line goes from 1990 to 1999. The corresponding correlation 
coefficients are shown in Table 13. 



mately 20%, in good agreement with the amplitude of the 
integrated measurement. Moreover, both contributions are 
in phase, which contrasts with the results found for the fre- 
quency shifts, where even modes seems to respond later 
than the odd ones. Odd and even TVP appear also bet- 
ter correlated around the maximum than during the low 
activity phases: both present exactly the same bump in 
the middle of the maximum of activity close to 1990, anti- 
correlated with decreasing activity during the same pe- 
riod, but their short term variations are different around 
1986 and 1996 during the activity minima 

Linear and rank correlation coefficients were calcu- 
lated with the same activity indicators than before and 
they are shown in Tabic The values found are large, 
showing a good correlation, but they are systematically 
lower than those for the frequency shift (range [0.67-0.84] 
against [0.86-0.94]). Moreover, the correlations are bigger 
for £=0,2 than for £=1,3 contrary to the results found for 
the frequency shifts. The linear correlation coefhcient be- 
tween frequency shift and TVP change (rp =-0.75, see also 
Table. ^ is well below the 0.9 reached in average between 
frequency shift and the solar activity indices suggesting 
that they are not linearly correlated and Fig. ^ shows in- 
deed that they tend to follow an hysteresis shape, rather 
than a strict line, when plotted against each other. 

This is important because if the decrease in TVP is due 
to the presence of local or surface activity, as we believe 
is the cause for the frequency shift, then they should be 
well correlated with almost no hysteresis. Therefore, the 
fact that the TVP is less well correlated with the surface 
activity indices and shows an hysteresis behavior when 
plotted against the frequency shift, may indicate that in- 
deed its variation is due to a decrease in the excitation 
efhciency or an increase of the damping rate at maximum 
which reflects a change in the convection zone structure 
that does not have to be correlated or strictly in phase 
with the surface magnetic features. The process of absorp- 
tion and damping of p-modes by an increasing number of 
rising flux tubes during the period of high activity ex- 
plored by Bogdan et al. (199^ ) is qualitatively compatible 
with our results. On the other hand, if the TVP varia- 
tions are only due to geometrical effects, this parameter 
would probably show a better correlation with the total 
area covered by the magnetic structures than with just the 
number of them. Although it is beyond the scope of this 
work, this should probably be further investigated to get 
a better picture of the possible sources of this phenomena. 
A more quantitative work including separate observations 
of the damping rate and amplitude variations of individ- 
ual modes ( Chaplin et al. 2000 ) is certainly needed to be 
more conclusive. These inferences from individual mode 
fits are potentially more informative but remains however 
less robust than those obtained from the fit of the cross- 
correlation of the power spectra. 
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Table 3. Correlation coefficients between the TVP and the solar indices where the TVP is defined as the product 
between the width and amplitude of the cross-correlation function between a spectrum and the 1986 reference spectrum. 
TVP refers to the average over the whole observed £ range [l — 0, 1,2,3) whereas TVPo,2 and TVPi.s refer to the 
average for even and odd degrees respectively, all quantities being frequency integrated values between 2.5 and 3.7 niHz. 
The last line gives the correlation coefficients between the TVP and frequency shift, is the Pearson linear correlation 
coefficient, rs the Spearman rank correlation coefficient and Pg is the probability of having no correlation. 



Index 




TVP 








TVPo,2 






TVPi.3 








rp 


rs 


Ps 




rp 


rs 


Ps 


rp 


rs 


Ps 




Ri 


-0.84 


-0.82 


4 10" 


8 


-0.87 


-0.88 


1 10-'' 


-0.79 


-0.71 


1 10" 


5 


Fw 


-0.83 


-0.79 


2 10" 


-8 


-0.86 


-0.87 


5 10-1" 


-0.78 


-0.69 


3 10" 


5 


KPMI 


-0.78 


-0.75 


2 10' 


-6 


-0.81 


-0.81 


7 10"** 


-0.74 


-0.67 


6 10" 


5 


MPSI 


-0.82 


-0.80 


3 10" 


-7 


-0.83 


-0.85 


9 10"" 


-0.80 


-0.69 


5 10" 


5 


TSI 


-0.76 


-0.67 


1 10" 


-4 


-0.80 


-0.76 


2 10"^ 


-0.70 


-0.58 


1 10" 


3 


He 


-0.84 


-0.82 


2 10" 


-8 


-0.87 


-0.87 


4 10-'° 


-0.77 


-0.72 


8 10" 


6 




-0.75 


-0.73 


4 10" 


-6 

















5. Frequency dependence of the frequency shifts 

In the previous sections, the variations of the frequency- 
integrated velocity power and frequency shifts have been 
studied. However, the time variation of the frequency 
shift is expected to be different at different frequencies. 
Previous works (e.g. Libbrecht fc Woodard 199"c| ; Anguera 
Gubau et al. 1992; Chaplin et al. 1998) have shown that 
the modes at high frequency become more sensitive to 
the solar cycle. The quality of our time-series and a new 
method to fit all the spectra at once motivates us to de- 
termine the frequency-dependence of the frequency shift 
for low degrees. 

The first method used is related again to the fact 
that pairs of low-degree acoustic modes are equally spaced 
in the spectrum. Each spectra is divided in regions of 
135 /uHz containing a set of modes £=0,1,2 and 3. Then, 
every region is cross-correlated with the corresponding re- 
gion of the reference spectrum (corresponding to the so- 
lar activity minimum of 1986), and the method explained 
above (Sect. ||) is used to calculate the frequency shift. 
Finally, the frequency shift of each region is fitted as a 
linear function of the integrated radio flux at 10.7 cm. 

In the second method proposed here, we try to fit to- 
gether all the spectra at once. We have established that 
the central frequencies of the solar acoustic modes vary 
during the solar cycle and there is a strong linear correla- 
tion with any of the solar indices. So, in order to improve 
the statistics we can fit all the spectra together introduc- 
ing the frequency shift as a new parameter. Here, as well 
as in the first method, the radio flux was chosen because, 
according to Table |l|, it leads to the best linear correla- 
tion coefficients for the central frequency variations but, 
as quoted before, this choice is not crucial as all indices 
present similar correlations. We emphasize that, while the 
flrst method is faster, the second provides us not only indi- 
vidual frequency shifts but also valuable mode parameters, 
i.e. mode resonant frequencies corrected for the solar-cycle 
effects. 



As the structure of the power spectrum is complicated 
by the presence of the 11.57/zHz sidebands, modes close in 
frequency must be fltted simultaneously in order to main- 
tain the stability in the fits. Therefore adjacent {n,i) and 
(n — 1, £-1-2) peaks are fltted together. The multiplet struc- 
ture induced by the rotation and the temporal sidebands 
for each mode of the pair are also included in the model. 
If we label by p the 60 /iHz wide part of the spectrum 
including the adjacent {n,£) and {n — 1, £ + 2) peaks and 

the mean frequency of the pair, the model for this part 
of the spectrum of each time-series i can be expressed as: 



k 

£+2 

E 



k=~l 
£+2 



E 



af^l^n£(r^/2)- 



- C™)' + (rp/2)^ 



"[™|^(n-l)(<?-H2)(rp/2)^ 



('^-^(„-i)(.+2)J' + (V2)^J 



Bp, (5) 



with: 



K^m = ^nl + 5v{Vp)i^t+2 ' {FIq ~ F^q) + m ■ S„£ + fc_D • fc,(6) 

where: 

— Vn£, is the central resonance frequency of the mode 
(n, I) at solar minimum; 

— 5v{vp)i^i+2, is the searched frequency shift per solar 
index unit, assumed to be the same for all the compo- 
nents of the pair; 

— Vp ~ (i^nf + i^,i-i,f+2)/2 is the center of the 60/xHz wide 
interval of the spectrum fltted; 

— Sn£, is the synodic rotational splitting for a multi- 
plet taken to be constant for all of them and equal 
to 400 nHz; 

— FIq, is the averaged radio flux for the given series [Fi^ 
being the value at the solar cycle minimum); 

— i?p, is a constant background level at the fltted fre- 
quency interval; 

— An£ , is the power at the resonance for sectoral compo- 
nents m=(. 
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— Q;[„j|i is a given ratio Anim/^ni constrained to take 
the theoretical value, for an instrument without spatial 
resolution, calculated by Christcnscn-Dalsgaard (1989) 
for the special case of resonant scattering observation 
using potassium line; 

— Fp, is the linewidth of the components of the multiplet 
assumed to be the same for all components of the two 
multiplets; and 

— /3|fe| , is the ratio of the power of the sidebands = 1) 
to the central peak (|fc| =0). These are directly related 
to the duty cycle and can be estimated empirically. The 
duty cycles of the yearly time-series remaining close to 
their average value of 25%, this allows us to fix the 
amplitude of the sidebands to /3|i|=0.5 for (3o=l. 

We assume that each m-component is well represented 
by a symmetric Lorentzian profile (even though Thiery 
et al. ( WOO) have found a slight asymmetry on them) and 
that the individual m-components are independent (see 



Foglizzo et al. 1998). Bp, A„i, Fp, Vni and Svne are the 
parameters to be fitted (vector a, hereafter). Notice that, 
in Sect. |^, it has been demonstrated that the TVP de- 
pends on the solar activity. This indicates that time varia- 
tions may also occur for amplitudes and linewidths. Such 
variations have indeed been found recently from BiSON 
( phaphn et al. 2000] ) for low degree and GONG ( ^omm 



of the linewidths and a decrease of the amplitudes leading 
to decrease of the TVP in agreement with our result. We 
therefore tried to parameterize amplitudes and linewidths 
as a function of time in the way used for frequencies, but 
the fits turned out to be very difficult and in many cases 
did not converge. 

At higher frequencies {i^ > 3500 /iHz) the peaks get 
wider, the width being greater than their rotational split- 
tings (Fp ^ Sni)- Moreover, the linewidths get so large 
that they become comparable to, or bigger than, the small 
frequency separations between the pair {£, n), {i+2, n~l). 
Therefore the fit is made at frequency intervals (labeled 
by q hereafter) of 165 /iHz centered at i^q and containing 
one pair of even modes (n, 0)(ri— 1, 2) labeled hereafter by 
p = 1 and one pair of odd modes {n,l)(n — 1,3) labeled 
hereafter by p = 2. Only one Lorentzian is fitted for even 
and another one for odd modes. Thus the model becomes: 



M;(z.,a)- 



2 



Apq{Tpq/2) 



ptt(^-^;.)' + (W2)^ 



B, 



9' 



(7) 



with: 



Vpq = iypq+Siy{iyq) ■ {F^ 



10 



(8) 



Here the parameters a to be fitted are: Bq, Apq, Fp^, 
and 5v(va). 



'pq 



As [Woodard (19841) points out, in the case of obser- 
vations without spatial resolution, the power spectrum 
of the solar p-mode oscillations is distributed around the 
mean Lorentzian profiles with a probability distribu- 
tion with two degrees of freedom. Consequently, the power 



spectrum appears as an erratic function where an abun- 
dance of frequency fine structure can be found. For this 
type of statistics the joint probability function associated 
to the observed power spectrum = {X^{vj)} cor- 
responding to the time-series i is given by: 



N 



1 



exp 



M*(i/j,a) 



(9) 



i=i ^ ^' ' 

where N is the number of frequency bins in the interval 
considered (i.e. the 60 /iHz wide interval for v < 3.5 mHz 
or the 165 /iHz wide interval for higher frequencies), and 
M* is given either by Eq. (^ or Eq. (0) depending also on 
the frequency domain. 

The likelihood for the 30 observed spectra is then given 
by the product of the individual likelihood function for 
each year. Thus, it can be written as: 



L{a)=l[f{X'^). 



(10) 



i=l 



We then look for the vector a that will maximize the like- 
lihood of the observed spectra according to our model. For 
numerical reasons one minimizes 5, defined as the nega- 
tive logarithm of the likelihood function. 



et al. ^000) for higher degrees, both showing an increase S{a 



30 N 

EE 

t=l J=l 



(11) 



To minimize this expression, we have used a modified 
Newton method (Press et al. 1992). The initial guesses 



for the parameters are important to avoid local minima. 
In that aspect, the frequency shift appears to be the more 
sensitive parameter and the initial guess was taken from 
the results shown Table || independently of the frequency 
range fitted. In the case of some parameters (amplitude, 
noise and linewidth), the natural logarithm of those have 
been fitted and not the parameters themselves. Doing this, 
S follows a normal distribution near the minimum and the 
covariance matrix for the vector a can be approximated by 
the inverse of the Hessian matrix found at the minimum 
of S. The uncertainties on each fitted parameter are there- 
fore taken as the square roots of the diagonal elements of 
the inverted Hessian matrix. 

The results are summarized in Fig. |[ In the top fig- 
ure, the triangles represent the frequency shift integrated 
in bands of 135 /iHz using the first method. The frequency 
dependence is quite clear, the frequency shift being close 
to zero at 2 mHz and then increasing progressively with 
the frequency to reach approximately 2 nHz/RF (where 
RF stands for radio fiux units, namely 10~^^J/s/m^/Hz) 
at the center of the p-mode and a maximum of 4 nHz/RF 
around 3.6 mHz, where it drops fast. The solid line denotes 
the inverse mode m ass extracted from the solar model of 
Morel et al. (1997 ). It has been averaged over the same 



regions in frequency and what we show represent the best 
fit to the data. The results corresponding to the simul- 
taneous fits are also shown in the same figure. The fre- 
quency shift for the even modes are represented by crosses 
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Fig. 5. The slope of the frequency shift in its assumed 
linear dependence with the radio flux is plotted here as a 
function of the frequency for low-degree p-modes. The re- 
sults for all modes (triangles) are obtained using the first 
method described in the text while the separate fits for 
odd and even degrees are obtained using a simultaneous 
fit (Eqs. (|[), (H)). The solid Hne is the best fit of the inverse 
mode mass calculated for the same set of modes. The bot- 
tom figure shows the results obtained at high frequency 
using Eq. (|^). The estimation of the frequency shifts cor- 



of Eqs. (0) and (|) gives the straightforward correspon- 
dance between b and the averaged frequency shift per ra- 
dio flux units. We note however that the two estimates are 
not equivalent due to the variation of the p-mode energy 
across the five minute band. The cross-correlation func- 
tion weight more the peak with higher amplitude located 
around ^ 2.9 — S.OmHz while the present integration gives 
the same weight to all the modes. 

At high frequency, the first points corroborate the fre- 
quency shift obtained using the first method, where a 
sharp downturn is found. Then, at higher frequency, a 
large fluctuation appears. A similar feature was found by 
Anguera Gubau et al. (1992 ) and also foun d for intermedi- 
ate degree by Libbrecht fc Woodard (199C ) in the analysis 
of BBSO data. These authors observed that the sensitivity 
of the mode frequency shif ts show similar sharp d ownturn 
located around 3.8 mHz. Chaplin et al. (1998 ) also re- 



ported (on their Fig. 4) a sharp downturn with negative 
frequency shift of about —10 nHz/RF around 4.3 mHz but 
they did not find the first downturn found at 3.75 mHz as 
in our analysis. The better sampling of the Mark-I data 
set and the fact that this first downturn was obtained with 
the two different methods we used for low and high fre- 
quencies make us confident in this result; moreover, this 
result is also found by Anguera Gubau et al. (1992 ). 

The frequency dependence of the frequency shift has 
been addressed from the theoretical point of view in differ- 
ent works. The model developed by Goldreich et al. (1991 ) 
suggests that a combination of an increase of the chromo- 
spheric temperature and a chromospheric resonance can 
be responsible for the sharp downturn at high frequency 
followed by an oscillation while the progressive increase in 
the five-minute band can be interpreted as an increase of 
the filling factor of the sma ll scale photospheric magnetic 
fields. On the other hand, jjain fc Roberts (1993i |1996|) 



argue that the presence of a magnetic field in the chromo- 
sphere and a combination of temperature and magnetic 
field strength variations could, qualitatively, explain the 
observed frequency dependence of the frequency shift at 
both low an d high frequencies. We n otice that both mod 
els, those of Goldreich et al. (1991 ) and Jain fc Robert! 



respon ding to low degree using the first method as well (1996|) present a wavelike behavior at high frequency qual- 



as the best fit of the inverse mode mass have been again 
plotted in the bottom figure as reference. RF stands for 
radio flux units, namely 10~^^ J/s/m^/Hz. 



itatively similar to the one found in our analysis. The par- 
ticular model of Goldreich et al. (1991) is even able to 



reproduce quantitatively this result, including the upturn 
around 4.5 mHz but no evidence has been found for the 



while black circles represent the odd ones. Both confirm 



the previous results showing a similar frequency depen- 
dence. These results are also in agreement with the anal- 
ysis carried out by Chaplin et al. (1998| ) on BiSON data 
for frequency below 3.6 mHz. The mean < 5v{i') > in 
the 2.5-3.7mHz range can be estimated by integrating the 
best fit of the inverse mode mass to our results divided by 
the length of the frequency range i.e. 1.2mHz. This leads 
to a value of 2.66nHz/RF compatible with the slope h re- 
ported in Tab. || for the linear dependance between the 
integrated frequency shift and radio flux. A comparison 



required chromospheric resonance (Woodard fc Libbrecht 
1991) and Kuhn (1998| ) point out that the change in the 



photospheric magnetic fleld strength needed in this model 
is much higher than the one obtained from recent infrared 
splitting observations of the quiet region fleld strengths or 
MDI magnetograms. Thus, Kuhn (1998| ) argues that pho- 
tospheric magnetic fluctuations arc unlikely to be respon- 
sible for the observed frequency shifts and proposes in- 
stead that turbulent pressure and mean solar atmosphere 
stratiflcation variations resulting from entropy perturba- 
tions through the solar cycle may be the dominant pro- 
cess affecting p-mode frequencies. In this model, the asso- 
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■ 1=0,2 

A e=l,3 




10 



100 



Fig. 6. The cross shows the normalized frequency shift for 
£=1-99 integrated between 2.5 and 3.7 mHz using LOWL 
data, and the solid line denotes the best fit of the inverse 
mode mass for those modes. We have also plotted the re- 
sults corresponding to very low degree using Mark-I data 
(Tab 1^): the dashed horizontal bar is for £=0,1,2,3; the 
triangle for £=1,3 is placed close to £=1 because this com- 
ponent gives the main contribution; and the square for 
£=0,2 has been arbitrarily placed close to £=2. 

dated temperature fluctuations may originate from near 
the base of the convection zone but in order to explore 
those possibilities and locate the different possible pertur- 
bations one needs to invert the even splitting coefficients 
(see Dzicmbowski ct al. 2000|) and track any fluctuation 
in the sound speed inversion using long term observations. 
We intend to extend our work in that direction using the 
observations of both low and intermediate degrees from 
LOWL and Mark-I. 

6. Analysis of the ^-dependence of the frequency 
shifts 



The LOWL instrument (Tomczyk et al. 1995), located 
in Mauna Loa, Hawai, is a Doppler imager based on a 
Potassium Magneto-Optical Filter, and it has been collect- 
ing solar observations for more than six years. With the 
installation of a similar experiment at the Observatorio 
del Teide, a new network called ECHO (Experiment for 
Coordinated Helioseismic Observations) intends to con- 
tinue the solar observations for a complete solar cycle 
(Tomczyk et al. 2000) with an increased duty cycle. 
Recently, six years of data have been re-analyzed through 
a new pipeline prod ucing mode parameter s for low and in- 
termediate degrees ( Jimenez-Reyes 2001 ). We have used 
the mode frequencies given by this analysis to compare 
the £- and frequency-dependence of the frequency shift at 
low and intermediate degrees. 

Figure ^ shows the normalized frequency shift using 
LOWL observations for £=1 up to 99. It has been per- 
formed using those modes between 2.5 and 3.7 mHz, as 




1.5 2.0 2.5 3.0 3,5 
Frequency (mHz) 



4.0 



Fig. 7. Normalized frequency dependence of LOWL data 
averaged for modes from £=1 up to 99 in intervals of 
150 ^Hz. The results for very low degree shown in Figure]^ 
are also plotted. The solid lines represent the best fit of 
the inverse mode mass to both data sets. 



we did for low degree. The inverse mode mass, which was 
calculated as well for each one of the modes fitted and 
then averaged in the same way, follows remarkably well 
the results. The figure shows as well the £-dependence of 
low degree p-mode frequency shifts for Mark-I found in 
Sect. ^. The general trend confirms the £-dependence of 
the frequency shift, the sensitivity at high degree (£=100) 
being almost twice to that of low-degree. We notice that 
the LOWL £=2 and Mark-I £ =0,2 are in good agree- 
ment and significantly lower then the inverse mode mass 
curve. However the LOWL data error bars in 5v for £ = 1 
and £ = 2 overlap significantly and the small value for 
£ = 0, 2 could either be due to particularly small 5v in 
£ = 2 but also in £ = 0. It will therefore be important 
to check this results with independent observations in the 
future. This is important because if a geometrical effect 
purely related to the integrated disk measurements affect 
the even modes measured by Mark-I, this should not be 
seen for the LOWL resolved measurements. If confirmed, 
this lower value of the frequency shift for £=2 may there- 
fore have a physical origin and be the signature of a local- 
ized perturbation. 

Figure |^ reproduces the frequency dependence for low 
degree shown in Fig. |^ together with the frequency shift 
averaged for £'s between 1 up to 99 in intervals of 150 /uHz 
obtained from LOWL data. Again, the best fit to the in- 
verse mode mass is shown. The ratio between the two 
slopes in the inverse mode mass fits is 0.74 and is nearly 
equal to the ratio between the mean mode mass calculated 
for both mode sets (0.75) showing that the £-dependence 
of the frequency shift is again well described by the £- 
dependence of the inverse mode mass. However, although 
the sensitivity is higher for higher degree modes (LOWL 
data), the fit is worse than for low degree (Mark-I data); 
moreover, in this later case, while roughly consistent with 
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the sizes of the error bars, the scatter seems to be orga- 
nized as an oscillation on top of the fitted line. The higher 
sensitivity of the frequency shift at high degree seems to 
be in contradiction with what Rcgulo et al. (1994) pointed 



out, but in agreement with recent analysis of Chaplin et al 
(1998|)] While all results point towards the existence of a 
perturbation confined close to the surface, there is still 
no convincing evidence of another cause that could exist 
deeper down as suggested by Regulo et al. (1994). On the 
other hand, the oscillation, also pointed out by Anguera 
Gubau et al. (1992), that seems to be present in the re- 
sults for low degree, has the same period 400-450 /iHz) 
that the one that is clearly found at higher frequencies 
(see Fig. || bottom) and deserves further attention and 
confirmation. 

Finally, we notice that the frequency dependence of 
the frequency shift obtained from the low degree modes 
of LOWL data is in very good agreement with the result 
plotted here using Mark-I data. This, added to the fact 
that the 6 years of LOWL data were covered also by Mark- 
I observations, make us very confident in the validity of 
our analysis of the i?-dependence using both instruments. 

7. Conclusion 

An analysis of the low-degree p-mode frequency shifts over 
more than a solar cycle has been carried out. Time depen- 
dence of the resonant frequencies and the total velocity 
power have been studied with several, old and new, meth- 
ods, yielding a correlation and an anti-correlation with 
the solar activity cycle respectively which confirms previ- 
ous results. Moreover, quantitative parameters have been 
introduced which yield the sensitivities of such variations 
with respect to a given activity index, in our case the ra- 
dio flux at 10.7 cm. The existence of a hysteresis behavior 
between both parameters seems to imply different direct 
causes for them. 

Also, the frequency shifts have been obtained as a func- 
tion of frequency and degree. These studies make use of 
data taken from two experiments (Mark-I and LOWL) 
which confirm each other results, when they coincide (at 
very low i modes), and complement their findings other- 
wise. Indeed, the main source for the variation of the fre- 
quencies with the cycle is located near the surface while 
a secondary, deeper rooted source, seems to be very weak 
if it exists at all. Interesting details such as the oscillation 
found for very low-^ p-modes and the surprisingly small 
frequency shift for i even modes need further investiga- 
tion. 

The results obtained and methods used in this work 
will also be useful in the analysis of very long time-series 
covering more than half the solar activity cycle, since both 
resonant frequencies and, to a lesser extent, amplitudes 
vary cyclically with time, providing wrong results if only 
data taken at parts of the cycle are used. In the following 
paper (second part of this work) we will use these methods 
to minimize the effects of solar activity on these param- 
eters in annual power spectra, allowing an average of all 



data thus improving the signal to noise ratio. Using this 
power spectra we will be able to give a new estimation of 
the solar rotational splitting, which combined with LOWL 
data, will be used to infer the solar rotation close to the 



core. 



Finally we note that the hysteresis between frequen- 
cies of odd and even modes or between mode frequencies 
and magnetic flux is not addressed in this work based on 
yearly time series. These studies are however of consider- 



able importance (see e.g. Moreno-Insertis & Solanki 2000) 



but would require shorter time series and better data (e.g. 
higher duty cycle) in order to reliably confirm previous 
analysis. This should probably be addressed in the future 
using long term observations of low degree modes from the 
ground based net works BiSON a nd IRIS-f -I- ( i.e. IRIS + 
LOWL +Mark-I, [Salabert (200l|) ). 
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